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Plasmons in two-dimensional systems with non-parabolic band, such as graphene, feature strong 
dependence on electron-electron interactions. We use a many-body approach to relate plasmon 
dispersion at long wavelengths to the Landau's Fermi-liquid interaction parameters and quasiparticle 
velocity. This general relation is illustrated using a model with N fermion species, giving a power- 
law dependence for plasmon frequency vs. carrier density at leading order in large N. An identical 
renormalization is predicted for magnet oplasmon resonance. We propose that gate-tunability of 
plasmons in graphene can be used to directly probe the effects of electron-electron interaction. 



Plasmonics has emerged recently as an active direc- 
tion in graphene research 1-5 . Surface plasmons in 2D 
electron systems are propagating charge density waves 
in which collective dynamics of clouds of charge is me- 
diated by electric field in 3D. 6 The dual matter-field na- 
ture of plasmons is a key ingredient for many interesting 
and important phenomena 7 ' 8 . Plasmons in graphene dis- 
play a range of potentially useful properties, such as low 
Ohmic losses, a high degree of field confinement, and 
gate-tunability, which make these excitations of great 
interest for optoelectronics 9 ' 10 . Gate-tunability of plas- 
mons in graphene was demonstrated recently 11,12 . 

The goal of this article is to relate the density depen- 
dence of plasmons with the interactions in the electron 
system. Plasmon dispersion depends on carrier density 
due to several effects. It takes simplest form in the limit 
of a weak electron-electron interaction, 1-4 



where Ep ~ n 1 / 2 is the Fermi energy of noninteracting 
massless Dirac particles, and n is carrier density. Here 
k is an effective dielectric constant of the substrate and 
a long- wavelength limit is assumed, g < pf. Interest- 
ingly, plasmon dispersion is also affected by interactions. 
Renormalization of the u> vs. q relation, Eq.(l), due to 
the effects of electron-electron interaction was predicted 
in Ref.13, where perturbation expansion in a weak fine 
structure parameter a = e 2 /hv was employed. The re- 
sults of Ref.13 point to an interesting possibility to di- 
rectly probe the effects of electron interactions by mea- 
suring plasmon dispersion relation. However, because of 
strong interactions in graphene, a ~ 2.5, the weak cou- 
pling approximation cannot serve as a good model, and 
a different approach is required. 

Acknowledging the difficulty of modelling the strong- 
coupling regime, it is beneficial to approach the problem 
at a somewhat more general angle. Rather than attempt- 
ing to make predictions based on a specific microscopic 
models, one can ask if a relation between the dispersion 
and some other fundamental characteristics of the sys- 
tem can be established. Below we point out that such 



a relation arises naturally from the Landau theory of 
Fermi liquids. 14 This theory affords a general, model- 
independent framework to describe systems of strongly 
interacting fermions at degeneracy. The effects of in- 
teractions are described by Landau parameters, repre- 
senting a "genetic code" of the Fermi liquid (FL). Their 
values can in principle be predicted from perturbation 
theory if interactions are weak. For systems with strong 
interactions, however, the most reliable way to obtain the 
Landau parameters is to use their relation with experi- 
mentally measurable quantities, such as compressibility, 
heat capacity, spin susceptibility and dispersion of col- 
lective excitations. 

Our many-body analysis upholds the conventional 
square- root dependence uj ~ q 1 / 2 . We show that all the 
effects of interactions are accumulated in the prefactor, 

oj 2 =YXq, Y = (l + F 1 )v, A=^^ (2) 

where pp is Fermi momentum, N — 4 is the number 
of spin/ valley flavors, and a long- wavelength limit is as- 
sumed, q <C pp. Here Fi is the Landau interaction har- 
monic with m = 1, and v is the Fermi velocity renormal- 
ized by interactions. The dielectric constant k may have 
a q dependence describing image charges arising due to 
conducting boundaries such as gates. In this case, a sim- 
ple model n{q) = \{ei + €2 coth(<7<f)) yields an acoustic 
plasmon dispersion to ~ q. 15 The quantity A in Eq.(2) has 
units of frequency and depends only on the fundamental 
constants and carrier density via pf. 

Magnetic field turns the gapless plasmon mode into a 
gapped mode. The magnetoplasmon dispersion relation 
obtained by adding Lorentz force to the Landau FL dy- 
namics takes the following form: 

CJ 2 B (q)=^(q)+Y 2 (eB/cp F ) 2 (3) 

where u>o(q) is plasmon dispersion at B = given by 
Eq.(2). The size of the gap at q — scales linearly with 
B, with the prefactor that varies inversely with pp ~ 
v}l 2 . Notably, the magnetoplasmon dependence on the 
interactions is described by the same combination Y = 
(1 + Fi)v as that appearing in Eq.(2). 



The quantity Y governs the interaction dependence of 
plasmon dispersion. Measuring it is a function of carrier 
density can shed light on the role of electron-electron 
interactions in the system. This behavior is in sharp de- 
parture with that for plasmons in two-dimensional sys- 
tems with parabolic band dispersion, where Galilean in- 
variance leads to an identity for Landau parameters, 
(1 + F%)v = «o j 14 where t>o = mpp is Fermi velocity 
of noninteracting particles at the same density. As a re- 
sult, the value (1 + F%)v is independent of interactions, 
leading to the "universal" long-wavelength plasmon dis- 
persion in the parabolic case, Wg(g) = ^f^" q at B = 
and ui 2 B (q) = ui 2 (q) + cu 2 for magnetoplasmon (where m 
is unrenormalized band mass and ui c = eB/mc is the 
cyclotron frequency). These dependences carry no infor- 
mation on the quantum effects and neither on the inter- 
actions. 

The situation is quite different in systems with non- 
parabolic dispersion, such as graphene. The density de- 
pendence in Y arises because the values v and Fi de- 
pend on the phase volume of the states in the conduc- 
tion band contributing to renormalization. The energies 
of these states span the domain from graphene band- 
width down to the Fermi energy, which is gate-tunable. 
To illustrate this point, we analyze the limit of a large 
number of spin/valley flavors, AT 3> 1. In this case, as we 
will see, Eq.(2) yields a power-law dependence on carrier 
density, giving 

u 2 oc n^/ 2 q (4) 

with the exponent value identical to that found from 
one-loop RG analysis for velocity renormalization, j3 — 
_8_ 16-18 p or a non _i n teracting system, f3 = 0. Crucially, 

the power law n^ 1-13 ^ 2 describes the dependence on car- 
rier density not only near charge neutrality but also for 
all accessible n values. 

Measurements of the density dependence of a plasmon 
resonance were reported in Rcf.f f . The observed depen- 
dence approximately follows the relation ui 2 oc q, with 
the prefactor exhibiting an approximately linear depen- 
dence on n 1 / 2 . However, the range of densities in which 
the dispersion was measured was not wide enough to pro- 
vide a clear distinction between (3 = and /? ^ 0. An 
attempt to experimentally determine the RG scaling ex- 
ponents from transport measurements was made recently 
in Ref.19, where a systematic variation of the period of 
quantum oscillations with carrier density was interpreted 
in terms of Fermi velocity renormalization, giving a value 
/3 = 0.5-0.55. This value is considerably larger than the 
one-loop RG result, /3 = -^-^ ~ 0.2, see below. This 
discrepancy is not yet understood. 

We note parenthetically that the interaction effects 
are not expected to vanish in graphene bilayer despite 
the parabolic character of its band dispersion. Elec- 
tronic states in graphene bilayer are Dirac-like rather 
than Schrodinger like, and hence do not admit Galilean 
transformation. For plasmons in this material we there- 
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FIG. 1: Resummed Feynman graphs for the polarization op- 
erator n(q, ui). The non-quasiparticle contribution IIo(q, uj) 
and the FL ladder ~}2 n>1 n„(q, ui) are shown. Only the con- 
tributions IIi (q, u>) and II2 (q, ui) contribute to the low-energy 
plasmon dispersion, see text. 



fore expect a behavior similar to that in materials with 
nonparabolic band, described by Eqs.(2),(3). 



I. MICROSCOPIC FERMI-LIQUID ANALYSIS 

The goal of this section is to relate plasmon dispersion 
with the standard quantities such as Landau FL inter- 
action parameters and renormalized velocity. The anal- 
ysis proceeds by standard steps via resumming the lad- 
der contributions to the dynamical polarization function 
which account for the quasiparticle dynamics in Landau's 
FL framework. This leads to a polarization response, 
H(q,uj) ~ q 2 /ui 2 , describing plasmon excitations in the 
low- frequency and long- wavelength domain, uj <C Ep, 
q < p F . 

Charge carriers in graphene single layer are described 
by the Hamiltonian for N = 4 species of massless Dirac 
particles. In second-quantized representation the Hamil- 
tonian reads 

H = ^p^oo-pV'p.i + ?4i-d, (5) 
p,i 

?4l-el=2 ^(q)V'p +q , 4 V'p'-q J V'p' J V'p, l ,(6) 
1,P,P',i,j 

where i,j = 1...N, vq ~ I0 6 cm/s is unrenormalized 
Fermi velocity, and V^(q) = 27re 2 /|q|K is the Coulomb in- 
teraction with the dielectric constant k describing screen- 
ing by the substrate. Here tjjp , is a two-component spinor 
with the components describing the wavefunction ampli- 
tude on two sublattices of the graphene crystal lattice. 
The amplitudes associated with the two sublattices are 
usually referred to as pseudospin up and down compo- 
nents, with the (pseudo)spin-l/2 Pauli matrices in Eq.(5) 
acting on (pscudo)spinors ipp.i- 

Plasmons are collective excitations of 2D electrons cou- 
pled by the electric field in 3D. They can be described 
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microscopically using the density correlation function 

K(d,u) = |^( Pq (i)p q (io))e^ ( *- to) (7) 

where p q (i) = En,,'/'p, ! (Wp+q,i( f ) are Fourier harmon- 
ics of the total electron density. The quantity K is ex- 
pressed in a standard fashion 14 through geometric series 
involving the polarization function II(q, uj) defined as the 
irreducible density-density correlator, 



n(q,w) 
e(q,w) 



e(q, W ) = l-V(q)n(q,u;) (8) 



Zeros of the dynamical screening function e(q, uj) give the 
poles of K, defining plasmon dispersion. To obtain the 
dispersion from the condition e(q, uS) = we need to gain 
information on II(q, uj) from a microscopic approach. In 
the long-wavelength limit, q <C pf, w <C Ep, the behav- 
ior of the quantity n(q, uj) is dominated by excitations 
near the Fermi surface, which can be described in the FL 
framework. 

The microscopic approach which forms the basis of the 
FL phenomenology involves several standard steps. At a 
first step, we isolate the contribution due to the quasi- 
particle pole of the electron Greens function G(x — x') — 
—i (ifj(x)ijj^ (%')) near the Fermi surface, 



(9) 



The first term is a regular part of the Greens function 
behaving as a smooth function near the Fermi level. The 
second term is a singular contribution describing quasi- 
particles, which dominates the FL dynamics: 



G (sing) (e,p) 



* e ~ £(p) + «7Sgne 



(10) 



Here Z is a quasiparticle residue, 7 is a quasiparticle 
decay rate, and £(p) = v(p—pp) is a quasiparticle energy 
dispersion, with v the renormalized velocity. 

This general relation can be specialized to the case of 
graphene as follows. The Green's function for electrons 
in graphene has a 2 x 2 matrix pseudospin structure. By 
projecting on the conduction and valence band, it can be 
represented as 



G(e,p)=G < (e,p)P < +G > (e,p)P > 



(11) 



where P 



>(<) = (1 ± ere p )/2 are projectors for the two 
bands (here e p is a unit vector in the direction of momen- 
tum p). The quasiparticle excitations with low energies, 
which govern the low-frequency and long-wavelength re- 
sponse, reside near the Fermi level. Without loss of gen- 
erality, we assume n-type doping, so that the Fermi level 
lies in the upper band, Ep > 0. In this case, excita- 
tions from the lower band do not appear explicitly in 
the FL theory and lead only to renormalization of vari- 
ous parameters such as the effective interactions and the 
quasiparticle velocity. The quasiparticle pole in Eq.(9) 



therefore arises only from the upper-band contribution 
G> whereas the lower-band contribution G< can be ab- 
sorbed into the regular part G^ reg ^ . Below the subscripts 
> and < will be omitted for brevity. 

The next step, which is key for understanding the role 
of low-energy excitations, is the analysis of the polar- 
ization function Tl(q, uj) at small uj and q. This is done 
by identifying the contributions due to pairs of Greens 
functions with proximal poles (the "dangerous" two- 
particle crosssections), 14 which we write symbolically as 
G (sing) G (sing) _ Z 2 ^^. One can represent U(q,uj) as 
a sum of terms with different numbers of such contribu- 
tions, 



nfa,w) = n («,w) + n 1 (g,w) + n 2 (g,w) + ... (12) 

Ili(q,u) ^T^GGT^, Tl 2 (q,u) = T U GGT"GGT U ... 

The corresponding graphs are shown in Fig.l. Here we 
introduced so-called quasiparticle-irreducible quantities: 
the renormalized scalar vertex and the two-particle 
scattering vertex T w (see Fig. 2 b,c). These quantities 
absorb all non-quasiparticle contributions in the upper 
band as well as the inter-band processes and the contri- 
bution of the states in the lower band. 

We recall that the quasiparticle-irreducible quantities 
are distinct from the conventional irreducible quantities 
defined as sums of Feynman graphs that cannot be split 
in two by removing two electron lines 14 . For example, the 
quasiparticle-irreducible vertex is obtained by sum- 
ming all kinds of graphs except the ones with dangerous 
cross-sections. The vertex (T^) can be obtained from 
the conventional irreducible vertex Tq (7o) by the resum- 
mation procedure pictured in Fig. 2, where the hatched 
blocks represent contributions due to pairs of Greens 
functions save for G( sing ) G^ sins \ 

To analyze the dependence on uj and q in the long- 
wavelength limit, we exercise caution and employ the 
quantities T w and T u taken at small but nonzero fre- 
quency and momentum values. We therefore adopt 
an approach similar to that used in Ref.20: our 
quasiparticle-irreducible quantities correspond to Lut- 
tinger's w-quantities which are taken at finite ui and q. 
They are distinct from the conventional w-quantities 14 
obtained in the limit u, q — > 0, (uj 3> vpq). This distinc- 
tion, however, turns out to be inessential: Luttinger's 
w-quantities reproduce the conventional w-quantites in 
the limit u, q — > 0, which can be taken in arbitrary order 
since dangerous cross-sections were left out of the defini- 
tion of T w and r u . 

Proceeding with the analysis, we note that the de- 
pendence on uj and q is very different for Ho(q,co) 
and H n >i(q,uj). We will first analyze the contribu- 
tion Ho(q, ui). This quantity does not contain dangerous 
cross-sections which can generate a nonanalytic behavior 
at small ui and q. Taking Tlo(q, uj) to be analytic, we can 
represent it as 



n (q,u;) =A(u})+B(cj) 



_q_ 

PF 



(13) 
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b) 




c) 





FIG. 2: Feynman graphs for the "quasiparticle-irreducible" 
quantities T", V". Bold lines represent a full Green's func- 
tion while the thin lines corresponds to a singular part c( smg ) . 
The vertex Fo represents the conventional irreducible vertex, 
whereas hatched blocks represent products of two Green's 
functions save for the contributions Q( sln -s) Q( sln -s) w hich give 
a non-analytical behavior at small w, q, see text. 



where A(ui) and B(oj) are regular functions. We recall 
that gauge invariance prohibits any physical response to 
spatially uniform time-dependent scalar field. Applying 
this to the full polarization function, Eq.(12), we con- 
clude that setting q — yields II (u>, 0) = 0. Also, since 
the contributions of the dangerous cross-sections GG van- 
ish at q = 0, all the quantities II„>i(g, u>) do so. We 
therefore conclude that the function A(uj) vanishes, leav- 
ing us with n (q, w) = {q/pp) 2 B{u}) + 0{q i ). This gives 
an effective q-dcpcndcnt permittivity 



k(q,u) = l-V(q)IL (q,u;). 



(14) 



Since V(q) = 2Txe 2 jnq in 2D, the effective long- 
wavelength permittivity is equal to unity and may hence- 
forth be ignored in the limit q/pF 0. 

It is instructive to compare the behavior of k(q, ui) with 
that arising for q pp. In this case the effects of finite 
doping are negligible and we can estimate the polariza- 
tion function using the result obtained for massless Dirac 
particles at zero doping, 



n(g,w) 



N 



16 y^2" 



(15) 



where N = 4 is the number of spin/ valley Savors. In 
the limit qv 3> u>, we obtain a well known q- independent 
result for permittivity 



k(q,u>) 



1 



irNa 



(16) 



This expression describes the effect of intraband polariza- 
tion in undoped graphene. We stress, however, that while 
the above expression is different from unity, it is obtained 



for q and uj which are not relevant for plasmon excita- 
tions. In fact, plasmons do not exist for such II(q, u>), 
as the plasmon dispersion terminates for q > pp. At the 
same time, the permittivity found in the long-wavelength 
limit g < pj? , w < , which is relevant for plasmons, 
equals unity. 

Next, we proceed with the analysis of the remaining 
terms, H n >\(q,u), which contribute to the low-energy 
plasmon dispersion. This is the part of polarization 
which account for the quaisparticle contributions. The 
corresponding Feynman graphs are given by a ladder with 
rungs consisting of two quasiparticle lines separated by 
vertex parts, as shown on a Fig.l. This gives geometric 
series that can be easily summed up: 



n(q, (J) = N j> t^T u 



vZ 2 



ui - qv(l + F) 



qvT" 



(17) 



where F is an integral operator 

Ff{6) = vZ 2 j fV(u>, q, 9, 00/(00. (18) 

Here v = pp/(2nh 2 v) is density of states per flavor, 9 is 
an angle between p and q. For zero external momentum 
q^O the kernel of the operator F depends only on the 
angle between p and p'. In fact, we will need only the 
quantity T"(0, 0, 9, 9') =T U (9- 9'). 

The scalar vertex on the Fermi surface with small 
external frequency and momentum Li, vpq ep can be 
decomposed as 



T"(w,q,0) = T +Tiuj + T2qcost 



(19) 



where 7o = Z^ 1 by virtue of Ward's identity 14 . The 
linear terms might have been important, if assessed by 
power counting. However these contributions drop out, 
because for external frequency w and momentum q, the 
expressions in question contain both T w (w,q, 9) and 
T w (-w, — q, 9). This leads to a cancellation of the terms 
linear in uj and q. 

In fact, only first two terms of these series are rele- 
vant for long- wavelength plasmons with vpq <C us <C £f- 
Anticipating the square root dependence for plasmon fre- 
quency vs. wavenumber, we expand in qvjui to obtain 



ni 



n 



2 — 



d9^,, „. uZ 2 vq cos 9 ^, , . 

— T>, q, ^ (-«, -q, 

27T u> — vq cos 9 

d9 (vqcos9) 2 v v 2 q 2 

27 J 2 = 2^"' 

d9d9',,, n ,vZ 2 vqcos9 
5 T (w,q,0) 

vZ 2 vq cos 9' 



= vZ 2 T$ j> 



(20) 



r>,q,0,0O 



UJ 



-T w (-a,,-q,0O 



v 2 Z% 2V X i f^r-(0 - 9') cos0cos0' 
us 2 J (2-k) z 

? 2 

w 2 q 2 2 

2 U! Z 



d8 
2^ 



r w (0) cos0. 



(21) 



5 



The terms Il„>3, expanded in qv/u>, yield contributions 
which are higher order in q . The same is true for con- 
tributions arising from expanding T u , T w in powers of q 
and uj (with the exception for potentially relevant linear 
terms 71, T2 in Eq.(19) which merely cancel out). These 
terms are therefore not essential in the long-wavelength 
limit. Combining all the above results for n and II n >i, 
we find the long-wavelength asymptotic behavior for the 
net polarization function: 



1 v 2 n 2 

n( q , w )= + 

Fi =vZ 2 I ^F"(0)cos0. 



(22) 
(23) 



The quantity F\ also represents an eigenvalue of the inte- 
gral operator F corresponding to eigenvectors cos 9 and 
sin6>. We can write Ft = F[cos0], which gives a Fourier 
harmonic of the operator kernel identical to Eq.(23). 

Plasmon dispersion can now be obtained from the re- 
lation 1 - V(q)U(q,u;) = 0, giving Eq.(2). The effects of 
interaction between fermions are "encapsulated" in the 
quantity Y which equals Fermi velocity v in the absence 
of interactions and is renormalized to a different value in 
an interacting system. 

We note a subtle difference between the quantities F m 
used in the FL literature 14 and those used here. The dis- 
tinction arises due to the long-range character of the 1/r 
interaction. In our case the density-density interaction 
F(6 — 0') accounts for the effects due to exchange corre- 
lation but not for the Hartree effects. The Hartree con- 
tribution is expressed through the 1/r interaction taken 
at the plasmon momentum q, corresponding to the Feyn- 
man graphs which can be disconnected by cutting a single 
interaction line. These contributions are incorporated in 
the dynamically screened interaction, Eq.(8), and hence 
not included in the definition of above. This difference 
manifests itself in the sign of the interaction parameters. 
In Fermi liquids with short-range interactions, the Lan- 
dau parameters describing density-density response are 
dominated by the Hartree effects. As a result, they have 
positive sign for weak repulsive interactions. In contrast, 
our F m are negative, since they are dominated by ex- 
change effects. In particular, we expect F% < 0. The 
negative sign, expected from this general argument, is 
also borne out by a microscopic analysis at weak cou- 
pling, see below. 

We also note an interesting analogy between the ap- 
proach developed in this section and the analysis of su- 
perconducting Fermi liquids by Larkin and Migdal, 21 and 
Leggett. 22 Refs.21,22 were concerned with Fermi-liquid 
renormalization of the quantities such as superfluid den- 
sity in a metal with BCS pairing. Their analysis focused 
on the current correlation function which determines the 
response of current to vector potential, and followed simi- 
lar steps as in the above discussion of II (q, u>). The renor- 
malization effects were expressed through a combination 
of FL parameters, featuring a cancellation for a system 
with a parabolic band. 



II. DENSITY DEPENDENCE FROM ONE-LOOP 

RG 

In this section we derive plasmon dispersion for a sim- 
ple model describing strongly interacting Dirac particles. 
This is done by employing the renormalization group 
analysis developed in Refs. 16-18 . We treat the Coulomb 
interaction which mediates two-body scattering by ac- 
counting for dynamical screening in the random-phase 
approximation (RPA), giving 



u =ZW 

e(q,w) 



e(q,w) = l-V(q)n(q } w). (24) 



Here the quantity e(q, uj) which describes dynamical 
screening is identical to that introduced in the above dis- 
cussion of the dynamical density correlator, Eq.(8). Here 
n(q, uj) is the polarization function 2 



II(q,a;)=iV J2 \ F *M 



/(ek, s ) - /(ek+q, s ') 



iuj 



£ k+q„ 



•»0' 



(25) 

with the band indices {s, s'} = ± and the coherence fac- 
tors |Fk,k'| 2 = |(k', s'\k, s)\ 2 describing overlaps of differ- 
ent pseudospin states. The polarization function is a sum 
of interband and intraband contributions, n = Hi + II2, 
described by s' ^ s and s' = s, respectively. The fac- 
tor e(q, uj) describes the effect of intrinsic screening in 
graphcnc arising due to both the interband and intra- 
band polarization. Below we give expressions that will 
be used in our analysis (for a comprehensive treatment 
of the quantity n(q, uj) we refer to Ref.2). For undoped 

graphene, only interband transitions contribute, giving 

2 

II 1 (q, uj) = — , — - q This dependence is key for 

V v 2 c\ 2 +uj 2 

RG analysis. 16-18 

The full RPA analysis of log-divergent corrections 
to Greens functions and vertices was performed in 
Refs. 16-18 . Below we show the results for one- loop RG 
calculation for large N, the results for N ~ 1 are qual- 
itatively similar, however the mathematical expressions 
are more combersome. The renormalization group flow 
for the velocity takes the form 



dv 



= pv, (3 = 



Nit 2 ' 



(26) 



where £ = ln(po/p) is the RG time parameter (here the 
UV cutoff is set by interatomic spacing in graphene lat- 
tice, po ~ a -1 ). This gives a power-law dependence 



«(P) = (p/Po) P v . 



(27) 



For TV = 4 we find (3 m 0.2. This value is obtained from 
a one-loop RG which employs 1/N as a small parameter. 
Acknowledging an approximate character of the scaling 
dimensions obtained from one- loop RG, we shall leave 
the exponent (3 unspecified in the analytic expressions. 

In the case of interest (doped graphene) the interband 
contribution III follows the above dependence for large 
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momenta and frequencies, |q| 3> pf, w 3> Ep, which 
provide dominant contribution to the RG flow of vari- 
ous quantities. The intraband contribution ±T 2 is much 
smaller than III at large energies and momenta, how- 
ever these two contributions are comparable for |q| ~ pp, 
uj ~ Ep. In the static limit, u> <C Ep, the polarization is 
dominated by the ±1 2 contribution. In the range q < 2pp, 
which is where we will need it below, it is identical to that 
for two-dimensional systems with parabolic band, 



n(|q| < 2p F ) = -Nv, 



(28) 



v = pp/2ttv (we refer to Ref.2 for the analysis of other 
regimes). This gives a standard expression for the static 
RPA-screened interaction 



2ne 2 



n\q\ + 2-nNve 1 



(29) 



We can obtain the two-particle scattering vertex I™ by 
taking the interaction on the Fermi surface, q ~ pp, 
ui Ep. This gives 

r»(9,9') = -.g P ,p'[r(Ap)] 2 f/ Api0 , Ap = 2 PF sin— , 

(30) 

where (7 PiP ' = |(p'a;'|pa:)| 2 = cos 2 (A6'/2) is the coherence 
factor describing the overlap of (pseudo)spinors describ- 
ing quasiparticles at different points of the Fermi surface 
(here A9 = 9 — 8'). The minus sign in Eq.(30) arises be- 
cause this expression represents a contribution from an 
exchange part of the two-particle vertex. 14 

Landau parameters can now be obtained from the re- 
lation between the Landau functional and the vertex T u , 
Eq.(18). Combining Eq.(18) and Eq.(30), we find 

F(9 - 9') = -iyg P ^Z 2 [T(Ap)} 2 U Ap . (31) 

In the large N limit, the static RPA-screened interaction 
can be approximated as J7 qjLJ ~ — n ^ ^ = , where we 
take into account that \Ap\ < 2pp. 

Both Z and T flow under RG, however their product 
remains equal to unity because of the Ward identity. As 
a result, the Landau function does not undergo power- 
law renormalization. Starting from Z(p)T(p) = 1, where 
both Z(p) and T(p) are given by power laws drawn from 
RG, we set p — pp. This gives 



(32) 



Up to a remnant dependence on pp which may arise in 
the agle dependence due to the ratio T(Ap)/T{pf), the 
function F does not flow under RG. 

While the function F(9 — 9') is essentially independent 
of doping, the velocity has a power-law dependence on 
doping, v oc Pp^. Combining these results, we find a 
power law dependence for plasmon dispersion, 



Ul 



A\q\, A 



Ne 2 ppv(p} 
2nh 2 



'(l + FOocp^, (33) 



where v vs. p dependence is given by Eq.(27). This re- 
sult is valid for plasmons with long wavelengths, q <C Pf- 
The predicted power-law dependence of plasmon on car- 
rier concentration holds for both large and small densi- 
ties, except very near the neutrality point where spatial 
inhomogeneity and thermal broadening start playing a 
role. Plasmons tend to stiffen due to RG-enhancement 
of velocity, and to soften due to the negative sign of F\ . 
However, since F± does not flow under RG, whereas ve- 
locity v does, the net effect of interactions leads to stiff- 
ening of plasmon dispersion. The predicted dependence 
A oc rS 1 ^^/ 2 can be used for extracting the exponent /3 
from measurement results. 



III. MAGNETOPLASMON IN A FERMI LIQUID 

Below we will treat plasmons using Landau's transport 
equations, presenting the derivation in a way that will 
help to make connection to the microscopic derivation in 
Sec. I. We will first deal with plasmons in the absence of 
magnetic field, then proceed to add a B field. Some of 
the relevant quantites, such the Landau FL interaction 
F(8 — 6') has already been introduced and analyzed, here 
we re-introduce and discuss them again in the framework 
of Landau FL transport equation. 

In Landau's scmiclassical framework, the fundamental 
effect dominating the Fermi-liquid behavior is forward 
scattering wherein the whole system of interacting parti- 
cles acts as a refractive medium which makes each quasi- 
particle energy a function of occupancies of other parti- 
cles. This is described by so-called Landau functional, 



Hp) 



d 2 P 



/(p,p')i5n(p / ,r) 



(34) 



where 8n(v, t) accounts for deviation of quasiparticle dis- 
tribution from equilibrium. 

Since deviation from equilibrium occurs in a narrow 
band of states near Fermi surface, it is convenient to 
write the Landau functional by setting |p| = |p'| = pp 
and parameterizing the Fermi surface by a unit vector 
n = p/|p|. Introducing the dimensionless Landau inter- 
action F(p,p') = ^/(p,p'), where v = pp/(2nh 2 v) is 
the density of states per flavor, we obtain 



/d9' 
— F(p,p')«Sn(p',r). 



(35) 



Here eo(p) = v(p — pf) is linearized quasiparticle en- 
ergy, the angle 9' describes orientation of p', and Sh(p) 
is obtained by integrating 5n(jp) along the Fermi sur- 
face normal. The expression (35) can be treated as a 
Hamiltonian of one quasiparticle moving in a selfconsis- 
tent field of other quasiparticles. The equation of mo- 
tions can then be obtained from Hamiltonian formalism 
via dt/n — {H, n}. This gives 

(d t + vV),5n(p, r, t) = -vVFfo(p, r, t), (36) 
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where F is the integral operator defined in Eq.(35). 

In a system with rotational symmetry, such as 
graphene and 2DEG's, the functional F depends only 
on the angle between p and p': 



F8n(9) 



2^ 



F(9 - 6')5n{6') 



(37) 



This expression defines a hermitian operator in the space 
of functions on the Fermi surface with the inner product 



</l(0)l/2(0)> = 



(16 
2^ 



(38) 



The eigenvalues of F are simply given by the Fourier 
coefficients 



F m = F{9)e- 



(10 
2^ 



F(6)e 



(39) 



The quantities F m represent the Landau interaction con- 
stants of a 2D FL system. 

To describe plasmons and relate their dispersion with 
the quantities F m we add to Eq.(36) a long range electric 
field arising due to oscillating charge density, 

(ft + vV(T+ F)) Sn(p,r,t) + eEV p n (p) = (40) 

where no(p) is the equilibrium Fermi distribution. Here 
E = — V$, where <£>(r) is the potential 

Here the sum is taken over N spin/valley flavors, and 
the dielectric constant k accounts for screening by sub- 
strate. Performing Fourier transform, Sn(9, r, t) = 

1 1 ^ty* 5 n tu,q(8) e ~ lu}t+t ' 11 ; we arrive at an eigenvalue 
equation of the form identical to that found in Sec. I by 
analyzing poles of the dynamical screening function, 



l-f(q)II(q,w)=0 



(41) 



The quantity n(q, w) is identical to that found above by 
summation of FL-type ladder graphs, 



1 



n(q,u) = NvTre - — ^qv 

U-qv(l + F) , 



(42) 



where trace is taken with respect to the inner product 
defined by Eq.(38). 

Plasmon dispersion in the long wavelength limit can 
be found by expanding in the ratio u|q|/o;. We obtain 



2 2 

vq v 



n(q,w) « -2 (qv|(l+F)|qv) = (cos 0\ (l+F)\ cosd) 

(43) 

Expressing the angle-averaged quantity through the 
Fourier coefficient F\ and using the relation v = 
Pp/(2ttv) we rewrite this result as 



H(4«) = ^«(l + *i) 



(44) 



Plugging this into Eq.(41) and restoring Planck's con- 
stant, we obtain the same expression for plasmon disper- 
sion as above, see Eq.(2). 

This analysis can be easily generalized to a system in 
the presence of an external magnetic field. This is done 
by accounting for the Lorentz force in the V p n term: 

ft + vV(T + F)j Sn(p, r, t)+(eE + -v x b) ■ V p n = 

where the velocity v = V p e(p, Sri) includes contributions 
due to the distribution function 8n(p). This equation can 
be linearized as above, rt(p,r,i) = no(p) + Sn(p,r,t). 
In doing so, particular caution must be taken with the 
Lorentz force term since it is affected by the FL inter- 
actions. Accounting for the term in the velocity v that 
depends on 5n(p), we write 

v(n + Sri) = d p e(n + Sri) (46) 
= <9 p e + dpFSn = v + Fdp'Sn, 

Here we used Eq.(35) and on the last step performed 

integration by parts. 

Terms linear in Sn can arise both from V p no and v. 

Taking a solution in a plane wave form Sn(r, p, t) oc 
e -iut+idT^ enQ ^f^ we have 

(iw - ivq cos 6>(1 + F)) f{9) - - (v x B) • V p /(0) (47) 

- (FVp/(0) x B) • v = ivV {q)qv cos 9 j> ^f(d') 

where 9 is the angle between q and v. This equation can 
be simplified as follows: 



iu - ivq cos 6(1 + F) - lj b (1 + F)d e j f(9) (48) 

■ vr \ a I d6 ' f(a>\ evB 
= ivV(q)qvcosO f — f(6 ), uj b = • 

J 27T p F C 

This is an eigenvalue problem with w a spectral parame- 
ter and f(9) an eigenfunction. Inverting the operator on 
the left hand side gives a self-consistency equation 



ivV(q)qv 



u + icon (l + F) d e - qvcos9(l + F) 



a dB i 

COS& — = 1, 

2tt 



Magnetoplasmon dispersion can be obtained via pertur- 
bation theory in the parameter qv/uj <C 1, giving 



W 2 (q) 



1 . 



(l + FO^ + a + FO-V^qV (49) 



Using the notation Y = (1 + Fi)v we arrive at Eq.(3). 
Magnetoplasmon dependence on interactions is described 
by the parameter Y identical to that found for plasmons 
at B = 0. The density dependence of the quantities v 
and 1 + F\ can be linked to their flow under RG. As 
discussed above, the power-law RG flow of velocity leads 
to stiffening of plasmon dispersion, which overwhelms the 
effect of softening due to the negative sign of F\ . 
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To put the above results in perspective, we recall that 
there is a crucial difference between the effects of interac- 
tion in graphene and in electron systems with parabolic 
band, such as those realized in GaAs. As noted above, 
for particles with quadratic dispersion the effects of in- 
teraction disappear because of the Fermi-liquid identity 
that relates renormalized velocity with the quantity Fi, 

Y = (1 + F{)v = v (50) 

where t>o = pp/m is Fermi velocity for noninteracting 
particles. This identity, which follows from Galilean in- 
variance, is peculiar to systems with quadratic dispersion 
and does not hold in graphene. As a result, plasmons in 
graphene depend in a nontrivial way on the Fermi ve- 



locity value v which undergoes renormalization by the 
interaction effects, and also on the FL interaction via a 
factor 1+Fi. These quantites, as we have seen above, are 
directly related to plasmon dispersion. The dependence 
of plasmon frequency on carrier concentration is expected 
to follow n 1 / 4 power law for weakly interacting electrons. 
The deviation from this dependence essentially maps out 
theoretically predicted RG flow of Y = (1 + Fi)v, with 
the momentum dependence parameterized by pp via the 
dependence on doping. This gives a unique opportunity 
to employ plasmon resonance as a probe of the RG theory 
of interaction effects in graphene. 

We thank F. H. L. Koppens, I. V. Kukushkin and M. 
Yu. Reizer for useful discussions. 
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